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Abstract 

The  development  of  wireless  communication  techniques  has  increased,  as  a consequence, 
the  use  of  positioning  systems.  Concerning  to  this,  GSM  providers  have  predicted  that 
position  location  will  be  an  integrated  service  in  the  3rd  generation  cellular  phones.  The 
techniques  that  are  used  by  each  one  of  the  current  positioning  systems  are  well-known; 
they  include  time-difference-of-arrival  measures,  angle-of-arrival  measures  and  power- 
of-arrival  measures.  In  this  paper  we  propose  a new  method,  which  is  based  on  the 
transmission  of  several  frames  with  different  frequencies  in  order  to  determine  the 
transmitter/receiver  distance. 


1.  Introduction 


Unlike  wired  communication  systems  that  are  static  and  easily  characterized,  radio  channels  are  highly 
random,  time  varying  and  difficult  to  predict.  The  propagation  path  in  a wireless  communication 
system  can  vary  in  real-time  from  close  range,  stationary,  line  of  sight  to  a heavily  obstructed  mobile 
path  due  to  both  moving  and  stationary  objects  [lj.  Needless  to  say,  although  wireless  communication 
systems  share  the  radio  channel  as  a common  limiting  factor,  the  channel  can  be  extremely  random  in 
nature  [2]  and  thus  plays  a major  role  in  the  design  of  many  radio  frequency  communication  systems. 

One  of  the  new  uses  of  these  systems  is  to  provide  accurate  user  position  location  (PL)  information, 
and  that  is  the  reason  why  several  different  position  location  technologies  have  been  developed  in  the 
last  few  years  [3],  As  it  is  well-known,  the  behavior  of  electromagnetic  waves  through  random  media 
is  modified  by  several  factors  that  appear  in  the  environment  [4],  If  we  know  the  limiting  factors  that 
restrict  the  propagation  of  the  waves,  the  medium  will  be  properly  characterized  and,  as  a result,  a 
relationship  between  the  propagation  velocity  and  the  frequency  of  transmission  will  be  defined. 

In  this  work,  we  propose  a new  method  to  determine  the  distance  between  the  transmitter  and  the 
receiver  by  solving  a system  of  non-linear  equations,  that  depends  on  the  physical  model  of  the 
medium  in  which  the  wave  is  propagated.  In  Section  2,  the  theoretical  basis  of  our  method  is 
described.  Consequently,  a system  with  a transmitter,  a microscopic  model  of  the  channel  and  a 
receiver  is  simulated  in  MATLAB.  The  results  are  shown  in  Section  3.  Finally,  some  conclusions  are 
remarked  in  Section  4. 


2.  Previous  Concepts  of  the  Multi-Frequency  Model 

The  propagation  of  electromagnetic  waves  depends  on  internal  factors  related  to  its  nature  and  external 
factors  that  are  characterized  by  means  of  the  medium  in  which  the  wave  is  propagated.  There  are 
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many  developments  in  this  field:  some  measures  of  the  GPS  propagation  principles  can  be  found  in 
[5],  where  reception  delays  due  to  atmospheric  phenomena  are  estimated.  Concerning  to  this,  the 
ionospheric  delays  are  considered  in  [6].  In  both  cases,  the  delays  not  only  depend  on  the  medium 
itself,  but  also  on  the  signal  frequency. 

Taking  into  account  these  studies,  a method  which  uses  a multi-frequency  technique,  that  is,  the 
transmission  of  several  frames  with  different  frequencies,  is  proposed.  Let  us  consider  the  general 
relationship: 

V/  = g(wi,kJ,...,kj)  (1) 


where  v,  is  the  propagation  velocity  and  depends  on  the  frequency,  g is  some  function  which  must  be 
determined,  w,  is  the  frequency,  and  ku...,kj  are  the  factors  related  to  the  medium  of  propagation. 
From  Eq.  (1),  the  distance  between  the  transmitter  and  the  receiver  can  be  estimated  by  using  a system 
of  non-linear  equations  [7].  These  equations  are  obtained  regarding  these  assumptions: 


• If  we  suppose  that  a frame  consisting  of  n frequencies  is  transmitted,  then  n equations  like  the 
following  ones  are  obtained: 


V|  'h  = d 


v„  • t„  = d 

n n 


(2) 


where  are  the  times  of  arrival  for  each  frequency  and  d is  the  distance  between  the 

transmitter  and  the  receiver. 


• If  the  time  difference  of  arrival  for  the  different  frequencies  is  calculated,  c„  then  we  will 
have  n-1  equations,  where  the  delay  between  each  one  of  the  frequencies  is  determined: 


ty  1 2 


t , —t„=c 


n-1 


(3) 


• Finally,  if  we  know  the  relationship  between  the  velocities/frequencies,  n equations  like  these 
ones  can  be  found: 

v,  = g(Wj  ,kx,...,kj) 


v„  = g(wn  ,ky,...,kj) 


(4) 


Therefore,  the  distance  d,  the  n propagation  velocities,  the  n times  and  the  factors  that  represent 
the  medium  must  be  calculated.  We  have,  then,  a system  with  3«-l  equations  and  2n+ 1 +j  unknown 
quantities,  and,  as  a consequence,  n=j  + 2. 


3.  Position  Location  using  the  Multi-Frequency  Model 

In  a common  way,  when  the  transmitter/receiver  distance  is  found  out,  the  equations  system  to 
determine  the  proper  position  of  an  object/person  can  be  obtained  as: 


{X,-Rj+(Y,-Ry}  + (z,-R,Y  = d,2 


(5) 
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where  (Rx,  Ry,  Rz)  is  the  receiver  position,  (Xh  Y„  Z,)  are  the  n transmitters  position,  and  d-,  are  the 
distances  to  each  one  of  the  n transmitters. 

The  parameters  of  the  model  that  characterizes  the  medium  can  vary  when  the  link  conditions 
change.  That  is  the  reason  why  the  positioning  system  will  have  to  update  the  model  in  real-time. 
Following  these  requirements,  a MATLAB  simulation  to  calculate  the  distance  in  a satellite-earth  link 
has  been  made.  Considering  the  Drude-Lorentz’s  model  for  the  ionosphere  [8],  we  have: 

.2 


k- 1 = 


2 2 
WQ  ~ W 


(6) 


where  kr  is  the  real  part  of  the  dielectric  constant  of  the  medium,  wp  is  the  plasma  frequency  and  w0  is 
the  natural  frequency.  If  the  refractive  index  is  n = , the  propagation  velocity  will  be: 


= g(w,kx,k2)  = - =■ 
n 


1 + - 


Wn 


■ w 


When  a two-parameter  model  is  used,  the  equations  system  will  be: 


(7) 


(8) 


where  i = 1, 2, 3, 4.  The  time  differences  can  be  calculated  as: 


h 


t2  - c, 


t]  t3  — c2 


U ~ C3 


(9) 


We  see  that  the  unknown  quantities  are  [f1;  t2,  t3,  t4,  ku  k2,  d).  As  a result,  we  can  deduce  that  the 
frames  must  be  sent  including  four  frequencies.  We  have  to  point  out  that  we  have  used  the  Newton’s 
method  to  solve  the  non-linear  equations  system.  That  is,  the  starting  equations  system  is: 


fl(xl,x2,x3,x4,x5)  = 0 


F(x)  = 


f7(xx,x2,x3,x4,x5)  = 0 


rfx(x)^ 


fl  (*) j 


= 0 


(10) 


Using  the  Taylor’s  expansion  of  Eq.  (10),  we  obtain: 


F(x)~  F(x0)  + J(x0)-(x-x0)  (11) 

where  J(x0 ) is  the  Jacobian  matrix  in  the  initial  point  x0=  [r10,  t2o,  t30,  t^,  k10,  k20,  dol 
As  a consequence,  the  following  conditions  must  be  considered: 

1.  A vector  with  some  initial  values  for  the  variables  of  the  system  must  be  introduced. 
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2.  After  that,  the  Jacobian  matrix  J(x 0(,_d)  for  the  current  solution  vector  3c0(M)  is 
calculated.  The  new  solution  vector  is  determined  as: 

*ot  = ^O(i-l)  — ^ (^00-1))'  ^(^O(r-l))  (12) 

where  J~](x00_}))  is  the  inverse  matrix  of  J(xou_ I}). 

3.  x0t  is  replaced  in  F(x)  and  the  error  value  is  obtained.  The  variation  of  the  error  is  revised 
to  establish  whether  the  solution  is  correct  or  not,  and,  then,  if  the  algorithm  is  finished. 

Table  1 shows  the  error  in  the  calculation  of  the  distance,  for  some  different  examples. 


Distance 

500 

1000 

5000 

Two-parameter  error 

1.2% 

0.1  % 

0.01  % 

Table  1:  Variation  of  the  error  and  distance 


4.  Conclusion 

In  this  paper  we  have  explained  a new  method  to  reach  a complete  position  location  service.  We  have 
shown  in  a theoretical  way  that  scattering  phenomena  can  be  used  to  estimate  the  transmitter/receiver 
distance  in  a radio  frequency  link.  Moreover,  any  electromagnetic  communication  in  which  scattering 
appears  could  be  considered  in  the  same  way.  Some  of  the  advantages  of  this  system  are,  on  the  one 
hand,  that  there  are  no  measures  in  absolute  time  and,  on  the  other  hand,  that  it  is  not  necessary  to 
synchronize  the  receiving  antennas  (see  [7]). 

As  a future  work,  we  propose  to  combine  this  method  with  a power-of-arrival  one  so  as  to  avoid 
the  wrong  distance  estimations.  Finally,  the  design  of  the  whole  positioning  system,  hardware  and 
software,  will  be  considered  as  well. 
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